Entanglement spectrum of a topological phase in one dimension 
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We show that the Haldane phase of S = 1 chains is characterized by a double degeneracy of the 
entanglement spectrum. The degeneracy is protected by a set of symmetries (either the dihedral 
group of TT-rotations about two orthogonal axes, time-reversal symmetry, or bond centered inversion 
symmetry), and cannot be lifted unless either a phase boundary to another, "topologically trivial", 
phase is crossed, or the symmetry is broken. More generally, these results offer a scheme to classify 
gapped phases of one dimensional systems. Physically, the degeneracy of the entanglement spectrum 
can be observed by adiabatically weakening a bond to zero, which leaves the two disconnected halves 
of the system in a finitely entangled state. 
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I. INTRODUCTION 

A topological phase is a phase of matter which cannot 
be characterized by a local order parameter, and thus 
falls beyond the Landau paradigm of condensed matter 
physics. Topological phases are typically characterized 
by a gap separating excitations from the ground state in 
the bulk and by the presence of gapless edge modes. The 
existence of edge excitations implies that a topological 
phase cannot be deformed continuously into a conven- 
tional, topologically trivial phase without going through 
a phase transition, in which the gap closes and the edge 
mode merges with the bulk. 

The Haldane phase of integer spin chainai*^ is an ex- 
ample of a "symmetry protected topological phase" in 
one dimension!^ This phase appears also in other one di- 
mensional systems, such as chains of interacting bosons^ 
and fcrmionsi^ These gapped phases lack a local order 
parameter, and arc not amenable to a description by a 
site factorizablc wave function. Alternatively, in certain 
cases, the Haldane phase can be characterized by the 
existence of fractionalized edge excitations, by a non- 
vanishing non-local "string" order parameter or, as re- 
cently proposed, by a quantized Berry phasci^^— 

However, in the most general case, the description of 
the Haldane phase in terms of a string order parameter or 
spin-i edge states is insufficient. As we will demonstrate 
below, slightly deforming the Hamiltonian can destroy 
the string order parameter (which is known to be fragile 
to small perturbations^) and lift the degeneracy of the 
edge states. Yet, as long as an appropriate set of sym- 
metries is preserved, the Haldane phase is stablo^ii^, in 
the sense that it is still separated from other, topologi- 
cally trivial phases by a thermodynamic phase transition 
in which the gap closes. This stability, by itself, can be 
used as an operational definition of the Haldane phase. 
However, it is desirable to find a definition which can 



be stated in terms of the ground state wavefunction of a 
single Hamiltonian. Recently, it has been proposed that 
topological phases can be characterized by their "entan- 
glement spectrum" , obtained by arbitrarily dividing the 
system into two parts, tracing out one half and diagonal- 
izing the reduced density matrix of the otheri^ir— This 
creates artificial edges, without disrupting inversion sym- 
metry. For example, the entanglement spectrum of the 
Affleck-Kennedy-Lieb-Tasaki (AKLT) stateii consists of 
two degenerate non-zero eigenvalues, which mimic the 
doubly degenerate energy edge spectrum of a system with 
a physical boundary.— 

In this paper, we show that the Haldane phase is char- 
acterized by a double degeneracy of the entire entangle- 
ment spectrum. This degeneracy is caused by the same 
set of symmetries which protect the stability of the Hal- 
dane phase, applied to the cigenstatcs of the reduced den- 
sity matrix. If the Hamiltonian is deformed while keep- 
ing these symmetries intact, the degeneracy remains until 
a phase boundary is crossed. This symmetry-protected 
double degeneracy can be used to define the Haldane 
phase in the most general situation, when both gapless 
edge states and a string order parameter are absent. 

The most surprising result of the analysis is that inver- 
sion symmetry alone is enough to preserve the degeneracy 
of the entanglement spectrum. If this is the only symme- 
try present, there are no gapless edge modes, since edges 
break inversion symmetry. There is also no string order, 
cither 4i 

This approach can be used to classify the phases of 
any one dimensional system, given its symmetry group. 
For a given set of symmetries, there are several types of 
gapped phases. One of them is the "non-degenerate" (or 
"trivial") phase, in which the eigenvalues of the density 
matrix can be non-degenerate. Besides this, there can be 
several types of "degenerate" ("non-trivial") phases. The 
entanglement spectrum in any one of the latter phases 
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has at least two-fold degeneracy. In this case, the density- 
matrix eigenstates transform in a non-trivial way under 
a projective representation of the symmetry group. 

The entanglement spectrum, although being associ- 
ated with a partition of the system at a certain point 
in space, actually carries highly non-local information 
about the ground state wave function. Wc show that 
the double degeneracy of the entanglement spectrum has 
a simple physical consequence. If one of the bonds of 
the system is adiabatically weakened until its strength 
reaches zero, the symmetry of the Hamiltonian continues 
to retain the degeneracy in the entanglement spectrum 
across the weakened bond. Hence, the von Neumann en- 
tropy of the partition in the final state is equal to ln(2) 
once the bond is broken. This is a physical reflection of 
the entanglement in the ground state, and can in princi- 
ple be used to identify it in experiment. 

This paper is organized as follows. First we introduce 
in Sec. |ll]a spin-1 model Hamiltonian which has a Hal- 
dane phase for a certain parameter range. In Sec. lIIIl we 
briefly review some properties of matrix-product states, 
which we use to study properties of the entanglement 
spectrum, and show how matrix-product states trans- 
form under symmetry operations. The main result of this 
paper, namely the degeneracy of the entanglement spec- 
trum in the Haldane phase, is derived in Sec IIVI Numer- 
ical results for several model Hamiltonians with different 
symmetries are shown in Sec. |Vl In Sec. IVII we briefly 
outline the generalization of these results towards a clas- 
sification scheme of gapped phases in one dimension. A 
more detailed discussion is deferred to a later publica- 
tion [isl A numerical experiment, which sheds light on 
the physical consequences of the degeneracy of the entan- 
glement spectrum, is presented in Sec. IVIII Finally, the 
results and conclusions are summarized in Sec. IVIIII 
Some details, concerning the application of the above 
results to generalizations of the Haldane phase and a 
derivation of the properties of the Haldane phase under 
inversion, are discussed in the appendices. 

II. MODEL HAMILTONIANS 

In order to study the stability of the Haldane phase, we 
will mainly focus on various spin-1 model Hamiltonians 
with different symmetries. As we will prove in Section lTVl 
the Haldane phase is protected by certain symmetries. 
One of them is a symmetry under a bond-centered spatial 
inversion 

(1) 

where Sj'^'^ are the spin-1 operators at site j. Other 
possible symmetries are the time reversal (TR) symmetry 

or the symmetry with respect to spin rotations by tt 
about a pair of orthogonal axes. As long as at least 



one of these symmetries is not broken, the entire en- 
tanglement spectrum remains doubly degenerate. There- 
fore, the Haldane phase maintains its identity and cannot 
evolve adiabatically to another phase. 

For concreteness, throughout most of this paper we 
consider the following spin-1 model Hamiltonian 

i j ] 

The symmetries of this model include translation, spatial 
inversion, a rotation by tt around the x-axis, and a com- 
bination of a rotation by tt around the y-axis and time- 
reversal e-^-^^xTR [which takes ^ Sf\ S^j 

-sy]. 

The phase diagram has been studied in Ref. [^. At large 
Uzz , we find a trivial insulator phase which can described 
by a caricature state where all the sites are in the \S^ =0) 
state. Furthermore, we find two antiferromagnetic phases 
Z2 and Z| with spontaneous non-zero expectation val- 
ues of (S*") and (5^), respectively. Uzz = = is 
the Hciscnberg point, for which one finds the gapped 
Haldane phase. Even when nonzero Uzz and values 
are introduced into this Hamiltonian, the Haldane phase 
is separated from the "non-degenerate phases" by phase 
transitions. Here, we investigate the question of what 
defines the Haldane phase and thus protects these tran- 
sitions from becoming smooth crossovers. We will later 
add different terms to Hq which break various symme- 
tries, and show explicitly for which classes of perturba- 
tions the Haldane phase remains well-defined. 

III. MATRIX PRODUCT STATE 
REPRESENTATION 

A. Definitions 

In order to prove the above statements, we will use 
a matrix product state (MPS) representation^^ of the 
ground state wavefunction. We also use this represen- 
tation to compute the ground state properties numeri- 
cally, using the infinite time-evolving block decimation 
(iTEBD) method,^ The iTEBD method is a descendent 
of the density matrix rcnormalization group (DMRG) 
method 1^ For the sake of completeness, we now review 
some of the properties of MPS's. A translationally invari- 
ant MPS for a chain of length L can formally be written 
in the following form^ 

l^-) = tr[r„,A...r„,^A] I mi...mL). (4) 

{rrij} 

Here, Tm are x ^ X matrices with x being the the di- 
mension of the matrices used in the MPS. The index 
m = —S, . . . , 5 is the "physical" index, e.g., enumerating 
the spin states on each site, and A is a x x Xj real, diagonal 
matrix. Ground states of one dimensional gapped sys- 
tems can be efficiently approximated by matrix-product 
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states22ri2ij in the sense that the value of x needed to ap- 
proximate the ground state wavefunction to a given accu- 
racy converges to a finite value as — > oo. We therefore 
think of X as being a finite (but arbitrarily large) number. 

The matrices F, A can be chosen such that they satisfy 
the canonical conditions for an infinite MPS ^^i^^ 



V r A^r^ = V A^r = i 



(5) 



These equations can be interpreted as stating that the 
transfer matrix 



(6) 



has a right eigenvector 5^^/ with eigenvalue A = 1. 
(* denotes complex conjugation.) Similarly, T^q,/.^^/ ~ 
Em(rm;3')*r^^A„A„, has a left eigenvector 8aa' with 
A = 1. We further require that ^^a' is the only eigenvec- 
tor with eigenvalue |A| > 1 (which is equivalent to the 
requirement that j'i/') is a pure state^i). 

The considerations given here become most intuitive 
when one considers, formally, an infinite chain. 

If the chain is infinite and has open ends, it may be 
partitioned at a certain bond. The wavefunction can then 
be Schmidt decomposed^^ in the form 



I*) =^A„|aL)|ai?), 



(7) 



where \aV) and |ai?) (a = 1, . . . , x) are orthonormal ba- 
sis vectors of the left and right partition, respectively. In 
the limit L — > oo, and under the canonical conditions 
([5]), the Schmidt eigenvalues Aq are simply the entries of 
the A matrix, K^a ■ are the eigenvalues of the reduced 
density matrix of either of the two partitions, and are 
referred to as the entanglement spectrum. The entangle- 
ment entropy is 5 = — J2a "^a ^a- This corresponds to 
the von Neumann entropy of the reduced density matrix. 
The states \aL) and \aR) can be obtained by multiply- 
ing together all the matrices to the left and right of the 
bond, e.g., if the broken bond is between sites and 



1, \aL) ^ E 



{m,},j<0 



nfe<0 ^^mk 



7a 



I . . . m^2m-imo) . 



Here, 7 is the index of the row of the matrix; when the 
chain is infinitely long, the value of 7 affects only an over- 
all factor in the wavefunction. Reviews of MPS's as well 
as the canonical form can be found in Refs. [26lf29l . 



u. 



(8) 



where C/s is a unitary matrix which commutes with the 
A matrices, and e*^^ is a phase. Thus, the matrices f/s 
form a dimensional (projective) representation of the 
symmetry group of the wavefunction. In close analogy to 
the derivation in Ref. [27l . we can derive a similar relation 
to Eq. ([8]) for time reversal and inversion symmetry. For 
a time reversal transformation Tm' is replaced by F^, 
(complex conjugate) on the left hand side. Finally, in 
the case of inversion symmetry F,„' is replaced by F^^, 
(transpose) on the left hand side of Eq. . 



IV. DEGENERACIES IN THE 
ENTANGLEMENT SPECTRUM 

We now turn to derive our main result, namely the 
degeneracies in the entanglement spectrum of the wave- 
function in the Haldanc phase. Our strategy is to de- 
termine when the transformation law for the Schmidt 
eigenstates under the symmetry operations of the system 
is non-trivial. From Eq. the Schmidt eigenstates of 
the left half of the system, \aL), transform under a sym- 
metry operation S according to the following rule: 



E|aL) =^(C/s)^a|/3L). 



(9) 



Similarly, the right Schmidt states \aR) transform by the 
conjugate matrix. Thus, the Schmidt eigenstates trans- 
form according to a projective representation of the sym- 
metry group of the system. The phases of the matrices 
Uyi arc not uniquely determined by Eq. ([S]), or by Eq. 

The phase ambiguities turn out to be the key to 
proving the degeneracies of the entanglement spectrum. 
We will show that for certain symmetries, there can be 
situations where the irreducible representations present 
in are all multi- dimensional. In these cases, which 
arc identified with the Haldane phase (or a generalization 
of it), the entire entanglement spectrum has non-trivial 
degeneracies. 



A. Inversion symmetry 



B. Symmetries in matrix product states 

In order to study the consequences of symmetries of 
the wavefunctions, it is useful to first study how these 
symmetries are reflected in the MPS representation. If 
1 4') is invariant under a local symmetry which is repre- 
sented in the spin basis as a unitary matrix then 
the F matrices can be shown to satisfy^! 



As a first example, let us consider a system which is 
symmetric under spatial inversion. The transformation 
law of F is written as 

rf„=e^'"t/]F,„C/j, (10) 

where Ux is a unitary matrix and 9x (z [0, 2tt) is a phase. 
Iterating this relation twice gives 
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T„,^e^''-{U^UiyT„,U^Ul (11) 
Now, the relation implies that 

TlAUjU^ATrn = e^'^'-UjUi, (12) 

m 

where we have used Eq. ((5)) and the fact that [Ui, A] — 0. 
Thus UjU^ is an eigenvector of the transfer matrix T [Eq. 
(|n])] with eigenvalue e^*^^. Since by our assumption, the 
only unimodular eigenvalue of T is A = 1 and this eigen- 
value is unique, we find that e^'^^ = 1 and Uj^Uj = e^'^^t 
where cj)! is a phase. Hence = C/ie~"^^. Repeating 
this relation twice, we find that e~^"^^ = 1, i.e. (pi = Q 

or TT. 

If 01 — TT, then Ui is an antisymmetric matrix. From 
this we find that all the eigenvalues Aq, are at least doubly 
degenerate. Moreover, the corresponding multiplicity ka 
is even for all a. This follows from the fact that Ui trans- 
forms the fcc-dimcnsional subspacc of states with eigen- 
value Aq, within itself. Therefore, the matrix (pro- 
jected into subspace a) satisfies detUj = det[(C/^)^] = 
det(-f/^) = (-1)''° detf/f . But since is unitary, 
detU^ ^ and therefore (-1)''° = 1. 

The fact that, in the presence of inversion symmetry, 
the phase 0x can only take discrete values (0 or tt) , leads 
to phase transitions between states when one would not 
expect them on the basis of the Landau paradigm of bro- 
ken symmetry. If an inversion-symmetric wavefunction 
evolves continuously, its characteristic phase 0x cannot 
change discontinuously, and therefore its value is fixed. 
The only way can change is through a critical point, 
where either the correlation length diverges because the 
transfer matrix T has a pair of unimodular eigenvectors 
(and the main relation, U'jUj = e^'^'t cannot be proven) 
or there is simply a discontinuous change in the ground 
state wavefunction (i.e., a first order transition). We can 
therefore identify two distinct states, characterized by 
= 0, TT. The state with (j)x = tt can be identified with 
the Haldane phase. To show this, we consider the AKLT 
state with = aa, A = "^^j where Ca (a = y, z) are 
Pauli matrices, and we use the time-reversal invariant 
spin basis \x) = ^ (|1) - |-1)), \y) = ^ (|1) + |-1)), 

\z) = |0). Under inversion, CTq — > crj — —cry(Ja<^y and 
thus Ui = Gy and 6'x = tt. Since cTj^c* = —1, we 

find that e"^^ = -1. The AKLT state is known to de- 
scribe the same phase as the Haldane phaseJi Therefore, 
we conclude that the Haldane phase is characterized by 
e'^^ = —1, e"^^ — —1, and a doubly degenerate en- 
tanglement spectrum. The wave function cannot evolve 
continuously if the phases Ox or 0x change discontinu- 
ously. This implies that changes of Ox or (px between 
and TT are always accompanied by a phase transition. 
Consequently, the degeneracy in the Haldane phase can 
only be lifted by a phase transition. The full argument 
for the existence of a transition in such case appears in 
Ref. S 



In the discussion above, we have assumed that the sys- 
tem is invariant under both inversion and translation [see 
Eq. (jl])]. However, in fact, inversion symmetry alone is 
sufficient to protect the double degeneracy in the entan- 
glement spectrum, as long as it is bond-centered. To 
show this, one can imagine adding a general commensu- 
rate perturbation to the Hamiltonian, such that the unit 
cell is enlarged. One can still write the ground state wave 
function in a translationally invariant form (j4]) were each 
site represents a single unit cell. If the unit cell is defined 
such that it ends at an inversion-symmetric bond, the 
new system is also inversion symmetric, and the entan- 
glement spectrum degeneracy remains protected. Since 
the size of the unit cell can be arbitrarily large, it is clear 
that translational symmetry cannot be essential for this 
argument to hold. The same argument can be made in 
the case of local symmetries, such as the ones described 
in Sections HVBIIV CI 



B. Time reversal symmetry 

The transformation of the MPS wavefunction F matri- 
ces under TR has the form 

E(^r)™™.r:„, =e'''-4F™[/^. (13) 

m' 

Here we have used the basis for the spins (to — 
— 1,0,1), and S7- = e^'^^ . From this, one can de- 
rive (in close analogy with the case of spatial inversion) 
that Uj-U^ = e"^^l where (pj- can be either or tt. 
If e*'^'^ = —1, then the double degeneracy of the en- 
tanglement spectrum follows (precisely as in the previ- 
ous section). The AKLT state F matrices transform as 
Pm — >■ cTyVm'^y. Thus Uj- = Gy and time reversal sym- 
metry is sufficient to protect the double degeneracy of 
the entanglement spectrum in the Haldane phase. 

The transformation of the MPS wavefunction F matri- 
ces under e^"^ xTR corresponds to a complex conju- 
gation (CC) of the wavefunction and has the form 

T*^ = e:'>^'^Ulcr..Ucc- (14) 

From this, one can derive that UceUcc = e'^'^'^t where 
(pec can be either or tt. The AKLT state F matrices 
transform as Fm — > Fm- It follows that Uqc — 1 and thus 
e'^'^^'xTR alone is not sufficient to protect the double 
degeneracy of the entanglement spectrum in the Haldane 
phase. Physically, this means that it is possible to add 
to the Hamiltonian a perturbation which is invariant un- 
der e"*'^'^ xTR but destroys the Haldane phase, in the 
sense that it is no longer separated from an unentangled 
product state by a phase transition. This perturbation 
has to break all other symmetries that may protect the 
Haldane phase (an example of such a perturbation can 
be found in Sec. |V|. 
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C. Sets of Rotations 

A symmetry of rotation about a single axis by 
where n is an integer, does not lead to any classification of 
phases. If E is a rotational symmetry of order n, one can 
show that J7£ = e**^ as in the previous section. Rescaling 

Uy: by e~ leaves Eq. (jSj) satisfied, and shows that ip has 
no significance. However, when there are multiple sym- 
metries, there is also a phase for each pair of symmetries 
Ei,E2. This phase is defined by noting that the trans- 
formation of Schmidt states corresponding to E1E2 can 
differ by a phase from the product of the Schmidt state 
representations of Si and E2: 

C/eiC/s. -e*''<''^'''=^C/siS,. (15) 

If the phase E2) cannot be gauged away by redefin- 
ing the phases of U^^ ^ > then the combined symmetry can 
lead to a protected Haldane phase. 

A concrete example is a system with symmetry under 
the dihedral group D2 of tt rotations about three orthog- 
onal axes (say, x,y and z axes). Since the product of tt 
rotations about the x and z axes TZ^T^z or TZzTZx give a 
TT rotation about the y axis TZy, the group is equivalent 
to Z2 X Z2.Thus it is sufficient to consider the action of 
two generators, say TZx and TZz- For TZ^i 

E = e^'-UlVr^U,, (16) 

where Tj^ = e"'^°'. Repeating this relation twice, we 
get Tm = e^**^- {UlfVmUl. From this it follows [analo- 
gously to the arguments below Eq. p^ ] that e^'^= = 1 
and = e**"=l. The phase factor e**'" is not impor- 
tant, since it can be absorbed in Ux- Therefore we 
can assume that Ul = 1. Similarly for U^, we arrive 
at = 1. The combined operation TZxTZz, however, 
may give rise to a non-trivial phase factor. By repeat- 
ing this symmetry twice, the associated unitary matrix 
UxUz can be shown (in the same way as above) to sat- 
isfy UxUz = e'^'^'^'UzUx- Since the phases of Ux and Uz 
have been defined, the phase factor e^'^'"' is not arbitrary, 
and can have a physical meaning. Clearly e'^'^'"' = ±1. If 
gi0a;j — then the spectrum of A is doubly degenerate, 
since A commutes with the two unitary matrices Ux, Uz 
which anti-commute among themselves. For the AKLT 
state, Ux = Ox and Uz = cr^, therefore UxUz — ~UzUx, 
and the Haldane phase is protected if the system is sym- 
metric under both TZx and TZz ■ 

V. EXAMPLES 

We now demonstrate how the symmetries discussed 
above stabilize the Haldane phase. We use the iTEBD 
method to numerically calculate the ground state of the 
model given by Eq. ([3]) , augmented by various symmetry- 
breaking perturbations. We used MPS's with a dimen- 



sion of X = 80 for the simulations. The double degener- 
acy of the entanglement spectrum is used to identify the 
Haldane phase. 
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FIG. 1: The colormaps show the entanglement entropy S for 
different spin-1 models: Panel (a) shows the data for Hamilto- 
nian Ho in Q , panel (b) for Ho plus a term which breaks the 
g-iirSf ^rpj^ symmetry [Eq. ^7}], and panel (c) for Ho plus 
a term which breaks the e"*''^ xTR and inversion symmetry 
[Eq. (|18p ]. The blue lines indicate a diverging entanglement 
entropy as a signature of a continuous phase transition. The 
phase diagrams contain four different phases: A trivial in- 
sulating phase (TRI) for large Uzz, two symmetry breaking 
antiferromagnetic phases Z| and Z|, and a Haldane phase 
(which is absent in the last panel). 

Example 1: We begin with the original Hamiltonian 
-ffo in (HI). This model is translation invariant, invari- 
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FIG. 2: Entanglement spectrum of Hamiltonian Ho in (|3]) 
for Bx = (only the lower part of the spectrum is shown). 
The dots show the multiplicity of the Schmidt values, which 
is even in the entire Haldane phase. 



ant under spatial inversion, under e 



and under 



xTR. Using the above argument, we know that 
inversion symmetry alone is sufficient to protect a Hal- 
dane phase. The phase diagram is shown in FIG. [T{a) 
and agrees with the results of Ref. A diverging en- 
tanglement entropy indicates a phase transition (see for 
example Ref. [3l|) and we use observables such as Sy, Sz 
to reveal the nature of the phases. In our approach, the 
Haldane phase can be easily identified by looking at the 
degeneracy of the entanglement spectrum as shown in 
FIG. [D the even degeneracy of the entanglement spec- 
trum occurs in the entire phase. 

In the entanglement spectrum of the TRI phase, there 
are both singly and multiply-degenerate levels. We have 
checked that the multiple degeneracies in the TRI spec- 
trum can be lifted by adding symmetry-breaking pertur- 
bations to the Hamiltonian (while preserving inversion 
symmetry). In the Haldane phase, on the other hand, 
the double degeneracy of the entire spectrum is robust 
to adding such perturbations. 

The colormap in FIG. |3{a) shows the difference of the 
two largest Schmidt values for the whole B^-Uzz phase 
diagram. In the Haldane phase, the whole spectrum is at 
least two-fold degenerate and thus the difference is zero. 

Example 2: The Hamiltonian Hq has in fact more 
symmetries than are needed to stabilize the Haldane 
phase. To demonstrate this, we add a perturbation Hi 
of the form 



qV qx' 
^3 . 



(17) 



11\ is translation invariant and symmetric under spatial 
inversion, but breaks the e"*'^'^ and the e^'"^' xTR sym- 
metry. The phase diagram for fixed Bz = O.IJ and 
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FIG. 3: The colormaps show the difference between the two 
largest Schmidt values |Ai — A2I for different spin-1 models. 
Panel (a) corresponds to the original Hamiltonian Hq in ([3|, 
panel (b) to plus a term that breaks the time reversal 
symmetry [Eq. (|17[l ]. and panel (c) to Hq plus a term which 
breaks time reversal and inversion symmetry [Eq. H18|l ]. The 
quantity |Ai — A2I is zero only in the Haldane phase. 



FIG 



= O.IJ as a function of B^ and Uzz is shown in 
[ijb). As predicted by the symmetry arguments 
above, we find a finite region of stability for the Hal- 
dane phase. This region is characterized, as before, by 
a twofold degeneracy in the entanglement spectrum, as 
shown in FIG.[2i;b). 

Example 3: In this example, we consider a case in 
which there is no symmetry that protects the Hal- 
dane phase. We add the following inversion symmetry- 
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breaking term: 

H, = 



3 

QZ ( QX QX 



■'j + l 

Note that this term is invariant under 



5!'5Li)+H.c.: 



(18) 



xTR. The 

phase diagram for the parameter R = O.IJ is shown in 
FIG. HKc). As predicted by the symmetry arguments 
above, we do not find a Haldane phase with a twofold 
degeneracy (see FIG. El^c)). The Haldane phase region 
is continuously connected to the TRI phase. The same 
scenario appears if we consider very small R. 

Example 4- Another example in which the Haldane 
phase and the TRI phase are continuously connected is 
recently given in Eqn. (6) of Ref. IH. Their model also 
does not have any of the symmetries which protects the 
Haldane phase. Thus their finding is consistent with our 
analysis. 



CLASSIFICATION OF GAPPED PHASES IN 
ONE DIMENSION 



VI. 



In Sec. IIVI we have identified several ^-parameters, 
such as 0x, (f'Tj Biiid (j)xz, which parametrize the phase 
ambiguities in the symmetry operations acting on the 
Schmidt eigenstates of the wavefunction. When one of 
these parameters is nonzero, the entanglement spectrum 
is degenerate, and a non-trivial (Haldane-like) "degen- 
erate" phase is stable over a finite range in parameter 
space. 

When more than one symmetry is present, the non- 
trivial ("degenerate") phases may be classified into sev- 
eral families depending on the combination of values 
taken by the corresponding 0's. In fact, there are even 
more phases than this argument would naively suggest 
Most generally, given the symmetry group of the system, 
the phases can be classified according to all the possible 
in-equivalent projective representations of the symmetry 
group. The general classification scheme of one dimen- 
sional gapped phases will not be described in detail in this 
work, but will be deferred to a later publication Ref. [l^. 

In Appendix. |X1 we consider various generalizations 
of the Haldane phase, which are protected by difi'erent 
symmetries. In a S* = 1 antiferromagnetic chain with 
Dzyaloshinskii-Moriya interactions in a magnetic field, 
we show that there is a stable Haldane-like phase which 
is protected by a modified inversion symmetry. The ex- 
tended Bose-Hubbard model also has a generalized Hal- 
dane phasci^ This phase is shown to be inequivalent to 
the usual Haldane phase of spin-1 chains, which is also 
supported by the symmetry group of this system. 



VII. ADIABATIC BOND WEAKENING 

The doubly degenerate entanglement spectrum is a 
unique feature of the Haldane phase, which can be used to 



0.8 

iog(2) 

0.6 
0.4 
0.2 




U =0.2J, B =0.3J 

ZZ X 

U =0.4J, B =0.3J 
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FIG. 4: Half-chain entanglement entropy of the model Hamil- 
tonian ((3)1 at a bond which is slowly weakened as a function 
of time Jwcak — J ~ tr and F — J^/40. The entanglement 
entropy of the resulting state is In 2 in the Haldane phase and 
zero otherwise. 



distinguish it from other phases. However, since the en- 
tanglement spectrum is a highly non-local property, one 
may wonder whether it has any physical consequences, 
which can be accessed in experiments. 

In the introduction, we discussed an adiabatic process 
in which a single bond in the system is slowly weakened 
to zero. The degree of correlation remaining across this 
bond can be measured by the entanglement entropy. We 
will show that, if the system starts in the Haldane phase 
and inversion symmetry about the weakened bond is pre- 
served throughout the process, the minimum value of the 
entanglement entropy of the two halves is ln(2). This is 
because the entanglement spectrum eigenvalues remain 
doubly degenerate. The minimum entropy is reached if 
just one pair of entanglement eigenvalues is nonzero. 

This property of the Haldane phase can be used, in 
principle, to identify it in an experiment. It means that, 
starting from the Haldane phase and separating it adia- 
batically into two halves A and B, some degree of en- 
tanglement between A and B must remain in the fi- 
nal state, as long as the symmetry that protects the 
Haldane phase is respected. This manifests itself in 
physical correlation functions between the two halves. 
Namely, there must exist a pair of physical operators 
Oa and Ob belonging to subsystems A and B, respec- 
tively, such that the disconnected correlation function 
Ca,b = {OaOb) — {Oa){Ob} remains non-zero, even 
when the two subsystems are completely disconnected. 
Starting from a "non-degenerate" phase, on the other 
hand, the final state can be completely unentangled 
across the cut bond.— 

To simulate the weakening of one bond numerically, 
we start by preparing ground states of Hamiltonian ^ 
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for different parameters using the iTEBD algorithm with 
a unit cell of size L = 80 which is large compared to 
the correlation length. We then evolve this state in time 
while decreasing the strength of one bond in between two 
half chains, Jwcak, according to Jwoak = J — tV. For the 
rate T = .P/iO and t = 0. . .40J"\ we found that this 
time evolution is essentially adiabatic. We calculate the 
entanglement entropy at the middle bond as a function 
of time. The result for Bx ~ 0.3 J and various values 
for Uzz is shown in FIG HI Within the Haldane phase 
{Uzz = 0.2J and U^z = 0.4J) the entanglement entropy 
at the end of the weakening process is equal to ln(2). 
In the TRl phase {U^z = 0.8J and Uzz = 1-OJ), the 
entanglement entropy decreases monotonically to zero. 

The robustness of the degeneracy in the entanglement 
spectrum of the Haldane phase has an intuitive explana- 
tion, as follows. Let us examine the Schmidt decompo- 
sition of the ground state wavefunction corresponding to 
dividing the system along the weakened bond. Initially, 
since = tt, the Schmidt states appear in doublets. As 
we show in Appendix B. the Schmidt decomposition can 
be written as 

x/2 

I*) = J2 ^2a-i(l«, - |a,2)|Za)), (19) 

Q = l 

where X2a-i are the Schmidt eigenvalues, \a,i) with i = 
1, 2 are Schmidt states of the left subsystem, and \a,i) = 
T\a,i) are the inversion-related states on the right sub- 
system. I^P) is odd under inversion, as can be seen by 
applying the inversion operator X. Since the Hamilto- 
nian remains symmetric under inversion throughout the 
bond weakening process, j^*) has to remain antisymmet- 
ric. Thus, at the end of the adiabatic evolution, the sys- 
tem is in the ground state of the antisymmetric sector, 
which generally differs from the true ground state (un- 
less the ground state of each of the disconnected halves 
is degenerate). This is because the true ground state 
is symmetric under inversion. The antisymmetric sector 
ground state can be written as --1=(|0)|1) — |1)|0)). Here 
|0) and |1) are, respectively, the ground state and first 
excited state of the left subsystem, and |0),|1) are the 
ground state and first excited state of the right subsys- 
tem, which are related by inversion to the corresponding 
states on the left. Thus the entanglement spectrum re- 
mains doubly degenerate, and entanglement entropy in 
the final state is In (2). 

Note that this property of the Haldane phase is not 
associated with the existence of zero energy edge states. 
(The two states |0) and |1) do not have to be degenerate.) 
In particular, the Hamiltonian Hq [Eq. ([3])] does not have 
any zero-energy edge state at an open boundary. 

VIII. SUMMARY 

In this work, we have considered the Haldane phase 
of S' = 1 chains as an example of a "topological" phase 



in one dimension. It has been known for a long time 
that this phase cannot be characterized by any local 
symmetry-breaking order parameter, and that its un- 
usual character only shows up in non-local properties, 
such as zero-energy fractionalized edge states and non- 
local string order parameters. When perturbing away 
from the SU(2) symmetric point, both the edge states 
and the string order can be eliminated. Remarkably, the 
Haldane phase can still remain stable, given that certain 
symmetries are preserved. I.e., the non-trivial topological 
character of the Haldane phase is protected by symme- 
try, even though the Haldane phase itself does not break 
any symmetry spontaneously. 

We have shown that the non-trivial "topological" na- 
ture of the Haldane phase of S" = 1 chains is reflected 
by a double degeneracy of the entire entanglement spec- 
trum. The degeneracy is protected by the same set of 
symmetries which protects the Haldane phase, and can- 
not be lifted unless either a phase boundary to another, 
"topologically trivial" phase is crossed, or the symmetry 
is broken. The Haldane phase is protected by any of the 
following symmetries: spatial inversion symmetry, time 
reversal symmetry or the dihedral symmetry D2 (rota- 
tions by TT about a pair of orthogonal axes). This result 
on the symmetry protection is completely consistent with 
what was obtained from different arguments.— The de- 
generacy of the entanglement spectrum can be used to 
characterize the Haldane phase in the most general sit- 
uation, in which edge modes and string order may be 
absent, (see TABLE 

The degeneracy of the entanglement spectrum in the 
Haldane phase is proven by examining how the Schmidt 
eigenstates transform under a projective representation 
of the symmetry group of the system. The transforma- 
tion laws contain phase factors, which are constrained to 
take discrete values by symmetry. If these phase factors 
are non-trivial, they require a degeneracy in the entangle- 
ment spectrum. Depending on which phase factors take 
non-trivial values, several distinct "Haldane-like" states 
arc possible. This offers a scheme to classify all possi- 
ble gapped phases of a one-dimensional system, given its 
symmetry group. Such a general classification will be the 
subject of a forthcoming paper 



symmetry 


string order 


edge states 


degeneracy 


D2 (=Z2XZ2) 


yes 


yes 


yes 


time reversal 


no 


yes 


yes 


inversion 


no 


no 


yes 



TABLE I: The different symmetries wliicli can stabilize the 
Haldane phase. For each class of symmetries, the table shows 
whether string order, edge states, or the degeneracy of the en- 
tanglement spectrum are necessarily present. The symmetry 
under tt rotations about a pair of orthogonal axes is repre- 
sented by the dihedral group D2. 

The degeneracy of the entanglement spectrum is a 



9 



highly non-local property, and is not easily related to 
physical observables. Nevertheless, the degeneracy of the 
entanglemerrt spectrum has direct physical consequences. 
It means that in the Haldane phase, the entanglement of 
the system across any cut cannot drop below the mini- 
mum value of ln(2). This can be observed, for example, 
by adiabatically weakening a bond to zero. In the Hal- 
dane phase, if inversion symmetry is preserved through- 
out this process, it leaves the two disconnected halves of 
the system in a finitely entangled state. In a "topolog- 
ically trivial" state, on the other hand, the two halves 
can be completely decoupled and form an unentangled 
product state after the process has ended. The non-zero 
residual entanglement is reflected in correlation functions 
of physical observables belonging to the two halves of the 
system. Such an adiabatic weakening process could, at 
least in principle, serve as a way to experimentally dis- 
tinguish the Haldane phase from other, "topologically 
trivial" phases. 



The following two cases often appear in models of mag- 



netism: a uniform DM interaction Z?, 



: D and a stag- 
Let us assume 



gered DM interaction Dj = (—lyu. 
D = (0, 0, D) (parallel to z-axis). 

The DM interaction, which is also known as anti- 
symmetric exchange interaction, clearly breaks inversion 
symmetry about the bond. Here we consider the Hamil- 
tonian of a 5' = 1 chain 



DM 



3 3 3 

(A2) 

where 77 = 1 for the uniform and ?y = — 1 for the staggered 
DM interaction case. This model breaks all the three 
symmetries we have discussed above, if i5 ^ 0. Thus 
one might expect that the Haldane phase is no longer well 
defined in this model. However, it turns out that the 
double degeneracy of the entanglement spectrum, and 
thus the well-defined Haldane phase, survives for D > 0. 
This can be simply understood because the Hamiltonian 



Acknowledgment 

We thank Ashvin Vishwanath, Ehud Altman, 
Emanuele Dalla Torre, Joel E. Moore, Michael Levin, 
Masaaki Nakamura and Yasuhiro Hatsugai for useful 
discussions. E. B. was supported by the NSF under 
grairt DMR-0757145. F. P. and A. M. T. acknowledge 
support from ARO grant W911NF-07-1-0576. M. O. 
was supported by KAKENHI grants No. 20654030 and 
20102008. The authors acknowledge the IPAM workshop 
QS2009 at which our work on the entanglement spectrum 
of the Haldane phase was initiated. 



Appendix A: Generalization of the 
symmetry-protected Haldane phase 

We have seen above that the Haldane phase of spin-1 
antiferromagnetic chains can be characterized by a dou- 
ble degeneracy of the the entanglement spectrum, which 
can be traced back to the non-trivial transformation law 
of its Schmidt eigenstates under certain symmetry oper- 
ations. The double degeneracy of the Haldane phase is 
protected either one of the three symmetries (dihedral 
group of TT-rotations about two orthogonal axes, time- 
reversal symmetry, or bond centered inversion symme- 
try). In fact, there are various generalizations of the 
Haldane phase, which are protected by modified sym- 
metries. 

For example, in magnetism, Dzyaloshinskii-Moriya 
(DM) interaction generally arises if the system lacks in- 
version symmetry about the center of bond. In the case 
of the one-dimensional chain, the DM interaction is given 

as 



(Al) 



H 



DM 



3 3 

B.Y.^3- (A3) 
i 

(A4) 



Here, Jl = ^/iP + D^ , if we choose 
for the uniform DM interaction, and 



(A5) 



(A6) 



for the staggered DM iirteraction, with a = tan~^ {D/ J). 
The resulting Hamiltonian is simply the standard XXZ 
antiferromagnetic chain in an magnetic field Bz- The in- 
version symmetry of the model guarantees a double de- 
generacy in the entanglement spectrum, and hence pro- 
tects the Haldane phase. 

In the context of the original Hamiltonian, however, 
the symmetry that protects the double degeneracy is 
somewhat obscured. The inversion T acts on the trans- 
formed Hamiltonian as 



H 



DM 



(A7) 



Here we define the irrversion I so that site j goes to site 
We find that the modified symmetry of the original 
Hamiltonian is the invariance under 



DM 



1^'^ HdmI' 



(A8) 



where 



1' = U(jIUl; = 



„iai:,(2j-l)S: 



(uniform DM) 
(staggered DM) 



(A9) 
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Namely, it is the invariance under inversion with an ap- 
propriate "twist" (rotation of each spin about z-axis). 

The invariance under X', which protects the Haldane 
phase, is not a generic symmetry and may be broken 
rather easily by perturbations which could occur natu- 
rally. For example, if the uniform magnetic field were 
applied to x-direction instead of z-direction in the Hamil- 
tonian (|A2|) . it is clear that the model is no longer invari- 
ant under X' . 

For a staggered DM interaction, a uniform field in x- 
direction leaves a staggered magnetic field cx 
after the transformation. The staggered field breaks in- 
version symmetry about the bond center and thus elim- 
inates the degeneracy in the entanglement spectrum. 
That a staggered field destroys the Haldane phase was 
noticed earlier 

As another example of physical interest, let us dis- 
cuss the "Haldane-Insulator" (HI) phase in the extended 
Bosc-Hubbard model (EBHM). We will show that the HI 
phase is protected by a similar mechanism. The model 
Hamiltonian of the EBHM reads 

Ubh = -f^(6]&j+i+H.c.) (AlO) 

3 

+ ^ XI " 1) + ^ X! 

3 3 

where we assume a filling of one bosons per site ((n) ~ 1) 
and t,U,V > 0. In Ref. [13, it has been shown that the 
EBHM has a phase which is analogous to the Haldane 
phase. This phase was termed a Haldane Insulator (HI). 

The symmetries of the EBHM are translation, time- 
reversal, inversion, and particle conservation. It is useful 
to consider an effective spin-1 model by truncating the 
Hilbert space of each site to states with n = 0, 1, 2, which 
is strictly justified in the large U limit. This modification 
is not expected to be important in this limit, since states 
with n > 2 are higher in energy. The corresponding 
effective pscudospin Hamiltonian reads 

3 3 

+Vj2S-S^+,+H' (All) 

3 

where we have introduced the pscudospin operator — 
n — I, and H' contains other terms which break the 
"particle-hole" symmetry of HaS, which is represented 
in the pscudospin language by a tt rotation about the 
X axis. This spurious symmetry is not crucial for the 
stability of the HI phase, as we shall show below. 

Ignoring H' , the Hamiltonian iJcff is very similar to Hq 
in Eq. ([3]), with the exception that the SJ'_^-^ + H.c. 
("hopping") term is of opposite sign. As a result, the HI 
phase of (jAll[) is not protected by inversion symmetry 
(X). The phase e*"^^, which is the parity of the ground 
state under inversion about a bond, is equal to -|-1 in 



this case; the ground state of a system with the ordinary, 
negative, sign for the kinetic energy cannot have nodes. 

However, the HI phase is protected by a modified sym- 
metry instead. The effective Hamiltonian can be mapped 
to the antiferromagnetic spin Hamiltonian ([3|) by a stag- 
gered rotation of spins by ±tt/2 about z axis, alternat- 
ingly on even and odd sites. The staggered rotation is 
given by the unitary transformation of the same form as 
eq. (|A6|) . but now with a = tt. We note that, if we in- 
crease the DM interaction from zero to infinity for the 
antiferromagnetic chain, \a\ changes from zero to 7r/2. 
Thus the present case is distinct from the antiferromag- 
netic chain with DM interactions. The transformation 
changes the sign of the hopping term to negative; in the 
spin chain context this makes in-plane exchange interac- 
tion antiferromagnetic as in Eq. ([3]). 

Following the discussion for a staggered DM interac- 
tion, and using a = tt, we find that the HI phase is 
protected by invariance under the operation 

As an interesting example, a staggered field in x-direction 
is now invariant under X' (and thus does not break the 
double degeneracy of the entanglement spectrum) , while 
a uniform field in the same direction is not. 

In terms of the discussion in Sec. IIVI , the HI phase is 
characterized by = tt and (pi = 0, and is thus distinct 
from the usual Haldane phase of Hq, with (px' = and 
(j>x = Ti"- This shows that these two states cannot be 
connected adiabatically while preserving either I or I'. 

Appendix B: Schmidt decomposition of a <j!>i = tt 
state 

Let us prove Eq. p^ . We consider an inversion- 
symmetric MPS 1^*) defined on a finite chain of an 
even length 2L and assume that j^P) is characterized by 
01 = TT. The MPS is written as 

I*) = E i4'^r„,A...r„,A 

{nij } 

xr„,^^jA. . .T„i^^VB}mi . ..m2L)- (Bl) 

Since we are interested in bulk properties in the limit L — > 
oo, we assume a sufficiently long chain with position- 
independent matrices r^. Vl and Vr are x dimensional 
column vectors which define the boundary conditions (to 
be specified later). Describing the boundary conditions 
in this way is possible as long as there are no edge modes, 
which is generically true when only inversion symmetry is 
present (Otherwise, the edges states of the two ends may 
require some extra care.). Since j^*) is invariant under 
inversion, the matrices r„i satisfy Eq. (jlOp . Applying 
this relation to the matrices Fm^^j^ . . . r„j2^, we get 

{mj} 
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xUiTl^^^A . . . Tl^,,uiVH\m, . . . m2L).(B2) 

Now, we choose boundary conditions such that Vr = 
UiVl- The wavefunction \^) can be written as 

\^)=e-^^0Y,K{UT)^p\am, (B3) 

where 

l")= E (^fr„,A...r„JJmi...mi), (B4) 

and \a) = T\a). Since (pi = tt, the Schmidt eigenvalues 
Act are all doubly degenerate (see Sec. II V A]) . Let us order 
the Ac's such that A2a-i = \2a for every 1 < a < %. The 
matrix Uap commutes with A. Therefore, it must have 



a block-diagonal form with 2x2 blocks on the diagonal. 
Since ~ —U (which follows from (/)x = tt, as shown in 
Sec. IIV Ap . the blocks on the diagonal of U are all of the 
form e*''°j(T2j where r]a is a phase. Therefore, we write 
l^") as 

x/2 

l*> = E 1)1", 2) - \a, 2)|a, 1)) , (B5) 

ct=l 

where \a,j) = e* 2 |2a — 1 + i) (j = 1,2) and 
\a,j) = T\a,j). In the limit L — > 00, the states \a,j) 
become orthonormal [as can be shown from the canoni- 
cal conditions ([S])], and therefore Eq. (|B5P is the Schmidt 
decomposition of \^). This concludes our proof. 
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